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Random sequential adsorption of straight rigid rods of length k (A;-mers) on a simple cubic lattice has been 
studied by numerical simulations and finite-size scaling analysis. The calculations were performed by using a 
new theoretical scheme, whose accuracy was verified by comparison with rigorous analytical data. The results, 
obtained for k ranging from 2 to 64, revealed that (i) in the case of dimers {k = 2), the jamming coverage is 
9j = 0.918388(16). Our estimate of differs significantly from the previously reported value of Oj = 0.799(2) 
[Y. Y. Tarasevich and V. A. Cherkasova, Eur. Phys. J. B 60, 97 (2007)]; (ii) Oj exhibits a decreasing function 
when it is plotted in terms of the k-mer size, being dj{o°) = 0.4045(19) the value of the limit coverage for large 
k’s; and (iii) the ratio between percolation threshold and jamming coverage shows a non-universal behavior, 
monotonically decreasing with increasing k. 

PACS numbers: 64.60.ah, 64.60.De, 68.35.Rh, 05.10.Ln 


INTRODUCTION 


The study of systems of hard non-spherical particles is one 
of the central problems in statistical mechanics and has been 
attracting a great deal of interest since long ago. In a sem¬ 
inal contribution to this subject, Onsager IjJ] predicted that 
rod-like particles interacting with only excluded-volume in¬ 
teraction can lead to long-range orientational (nematic) order. 
This nematic phase, characterized by a big domain of parallel 
molecules, is separated from an isotropic state by a phase tran¬ 
sition occurring at a finite critical density. Numerous model 
systems of anisotropic objects, ranging from inorganic to or¬ 
ganic and also biological particles, were analyzed by using the 
Onsager’s theory. An extensive overview of this field can be 
found in the excellent reviews of Refs. Hi. 

The long-range orientational order disappears in the case 
of irreversible adsorption (no desorption), where the distri¬ 
bution of adsorbed objects is different from that obtained at 
equilibrium. Thus, at high coverage, the equilibrium state 
corresponds to a nematic phase with long-range correlations, 
whereas the final state generated by irreversible adsorption is 
a disordered state (known as jamming state), in which no more 
objects can be deposited due to the absence of free space of 
appropriate size and shape. The jamming state has infinite 
memo^ of the process and the orientational order is purely 
local 

The phenomenon of jamming plays an important role in nu¬ 
merous systems where the deposition of objects is irreversible 
over time scales of physical interest 101. The random sequen¬ 
tial adsorption (RSA), introduced by Feder 101, has served as a 
paradigm for modeling irreversible deposition processes. The 
main features of the RSA model are: (1) the objects are put 
on randomly chosen sites, (2) the adsorption is irreversible 
and (3) at any time only one object is being adsorbed, so 
that the process takes place sequentially. In this theoretical 
frame, anisotropic particles of different shapes and sizes have 


been studied: linear k-mers (particles occupying k adjacent 
lattice sites) flexible polymers HI ll^. T-shaped ob¬ 

jects and crosses 


, squares Ill7l420ll . disks 112111 , regular and 
star polygons 1I22|] . etc. In all cases, the limiting or jamming 
coverage, 9j, strongly depends on the shape and size of the 
depositing particles. 

In the case of lattice models straight rigid rods (or linear k- 
mers), which is the topic of this paper, the RSA problem has 
been exactly solved for one-dimensional (ID) substrates. In 
this limit, 0(f) can be written as |0], 


0(f) = k f exp 
Jo 


f I 


du. (1) 


From Eq. (ID, the dependence on k of the jamming cover¬ 
age can be obtained. Note that 0(f) represents the fraction 
of lattice sites covered at time f by the deposited objects and, 
consequently, 9{t = °o) = 9j. 

In two and three dimensions, the inherent complexity of 
the system still represents a major difficulty to the develop¬ 
ment of accurate analytical solutions, and computer simula¬ 
tions appear as a very important tool for studying this subject. 
In this direction. Bonnier et al. investigated the deposition of 
linear k-mers (with k ranging between 2 and 512 on a two- 
dimensional (2D) square lattice H. The authors found that 
the jamming concentration monotonically decreases and tends 
to 0.660(2) as the length of the rods increases. 

Vandewalle et al. studied the percolation and jamming phe¬ 
nomena for straight rigid rods of size k on 2D square lattices 
H and found that, for values of k between 2 and 10, the 
ratio 9p/9j (being 9p the percolation threshold) remains con¬ 
stant 9p/9j 0.62, regardless of the length of the particle. 

Based on this finding, the authors suggested that both critical 
phenomena (percolation and jamming) are intimately related. 
Kondrat and Pfkalski |12] extended the study of Ref. H 
to larger lattices (lattice size L = 30,100,300,1000 and 2500) 
and longer objects (1 < k < 2000). The results obtained re- 
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vealed that: (1) as reported in Ref. the j amming coverage 
decreases monotonically approaching the asymptotic value of 
9j = 0.66(1) for large values of t, (2) the percolation thresh¬ 
old is a nonmonotonic function of the size k, showing a min¬ 
imum around k — 13; and (3) the ratio of the two thresholds 
0p/6j is nonmonotonic too: after initial growth, it stabilizes 
between k — 3 and k = 7, and then it grows logarithmically. 


Recently, Lebovka et al. studied an anisotropic RSA of 
straight rigid rods on 2D square lattices lll3n . In this model, 
the vertical and horizontal orientations occur with different 
probabilities, and the degree of anisotropy of the system can 
be characterized by an order parameter measuring the differ¬ 
ence (normalized) between the number of line segments ori¬ 
ented in the vertical direction and the number of line segments 
oriented in the horizontal direction. The authors investigated 
the effect of k-mer alignment on the jamming properties. In 
the case of isotropic systems (order parameter equal to zero), 
the results obtained by Lebovka et al. are in excellent agree¬ 
ment with previous simulation data in Ref. lll2ll . 


The RSA problem becomes more difficult to solve when the 
objects are deposited on three-dimensional (3D) lattices, and 
only very moderate progress has been reported so far. In the 
line of present work, Tarasevich and Cherkasova examined the 
percolation and jamming properties of dimers (straight rigid 
rods with k = 2) on 3D simple cubic lattices and found that 
(1) 9p = 0.2555(1) and 9j = 0.799(2), and (2) in contrast with 
percolation and jamming on square lattices, 9p/9j 0.32 for 

cubic lattices lE^ . 


In the present paper, the study of Tarasevich and 
Cherkasova is extended to larger k-mers {2 < k < 64). Us¬ 
ing computational simulations and a new method introduced 
here, the jamming properties of straight rigid rods on 3D sim¬ 
ple cubic lattices are studied. The method is based on the 
calculation of the probability Wl{9) that a lattice composed 
of Lx Lx L elements reaches a coverage 9. As it will be dis¬ 
cussed in details in the next section, the intersection point of 
the curves Wl{9) for different values of L gives an accurate es¬ 
timation of the jamming threshold in the infinite system. The 
main objectives of the paper are (1) to evaluate the accuracy 
and applicability of the new technique to estimate jamming 
properties, (2) to determine the dependence of the jamming 
coverage on the size of the deposited k-mers, and (3) to com¬ 
pare our results with the corresponding ones in Ref. 1231] . The 
work is a natural extension of our previous research in the area 
of percolation and jamming of linear k-mers. 


The outline of the paper is as follows. In Sec. II we describe 
the deposition model, along with the simulation scheme and 
the procedure used to calculate the jamming coverage. In Sec. 
Ill we present the results and discussion of the simulations. 
Finally, the general conclusions are given in Sec. IV. 


MODEL, BASIC DEFINITIONS AND SIMULATION 
SCHEME 


The following scheme is usually called standard model of 
deposition or RSA. Let us consider the substrate is represented 
by a three-dimensional simple cubic lattice of M = LxLxL 
sites and periodic boundary conditions. In the filling process, 
straight rigid rods of length k (with k > 2) are deposited ran¬ 
domly, sequentially, and irreversibly on an initially empty lat¬ 
tice. The procedure of deposition is as follows: (/) one of the 
three possible lattice directions {x,y,z) and a starting site are 
randomly chosen; (ii) if, beginning at the chosen site, there 
are k consecutive empty sites, then a k-mer is deposited on 
those sites. Otherwise, the attempt is rejected. When N rods 
are deposited, the concentration is 0 = kN/M. 

Due to the blocking of the lattice by the already randomly 
deposited dimers, the limiting or jamming coverage, 9j = 
9{t = oo) is less than that corresponding to the close pack¬ 
ing (0; < 1). Note that 0(f) represents the fraction of lattice 
sites covered at time t by the deposited objects. Consequently, 
0 ranges from 0 to 9j for objects occupying more than one 
site, and the jamming coverage depends on the size of the de¬ 
posited object fl&Sl- 

It is well known that it is a quite difficult matter to analyti¬ 
cally determine the value of the jamming coverage for a given 
lattice. For some special types of lattices, geometrical consid¬ 
erations enable to derive their jamming thresholds exactly i^. 
For systems which do not present such a topological advan¬ 
tage, jamming properties have to be estimated numerically by 
means of computer simulations. 

In order to calculate the jamming thresholds, we introduce 
the probability Wl{9) that a lattice composed of Lx Lx L ele¬ 
ments reaches a coverage 0. In the simulations, the procedure 
to determine Wl{9) consists of the following steps: (a) the 
construction of the M = Lx Lx L lattice (initially empty) and 
(b) the deposition of particles on the lattice up to the jamming 
limit 9j. In the late step, the quantity mj{9) is calculated as 


m{9) 


1 for 0 < 9j 
0 for 0 > 9j. 


( 2 ) 


n runs of such two steps (a)-(b) are carried out for obtaining 
the number m(0) of them for which a lattice reaches a cover¬ 
age 0, 


fn{9) = Y^mi{9). 

i=i 


(3) 


Then, Wl{9) = m{9)/n is defined and the procedure is re¬ 
peated for different values of lattice sizes. For infinite systems 
(L —oo), Wl{9) is a step function, being 1 for 0 < 9j and 0 
for 9 > 9j. For finite values of L, Wl{9) varies continuously 
between 1 and 0, with a sharp fall around 9j. In conclusion, 
Wl{9) appears as a good parameter evidencing the jamming 
threshold and its definition is computationally convenient. 

Using a system of fully aligned rods deposited on 3D cubic 
lattices as an example of application 0251], we show in the fol- 
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FIG. 1: Curves of Wj, vs 6 for fully aligned dimers on simple cubic 
lattices of different sizes as indicated. From their intersections one 
obtained 6j. In this case, 6j = 0.864656(18). 

lowing how to determine the corresponding jamming cover¬ 
age from the probabilities V^(0). This model of fully aligned 
rods can be exactly solved jot], which allows us not only to il¬ 
lustrate the use of the method introduced here, but also to test 
its validity and accuracy by comparing our simulation results 
with those of Ref. |@] . Five k-me.r sizes were chosen for the 
comparison k{= 2,3,4,5,6). Then, a set of n = 10^ runs were 
carried out for each pair k and L/k{= 18,24,30,36,42). 

Table I; Simulation and theoretical values of Oj (as indicated 
in the text) for fully aligned rods on a simple cubic lattice and 
k ranging from 2 to 6. 



Aligned rods 

k 

9j [Eq.(l)] 

Oj [Hi's] 

2 

0.864665 

0.864656(18) 

3 

0.823653 

0.823646(16) 

4 

0.803893 

0.803892(12) 

5 

0.792276 

0.792284(15) 

6 

0.784630 

0.784637(16) 


In Fig. 1, the probabilities Wl( 0) are presented for a sys¬ 
tem of aligned dimers {k = 2) and several values of L as indi¬ 
cated. As it can be observed from the figure, curves for dif¬ 
ferent lattice sizes cross each other in a unique point, which is 


FIG. 2: Curves of WT as a function of the density 6 for two values of 
fc-mer size (k = 2, curves on the right, and k = 5, curves on the left) 
and lattice sizes ranging between L/k= 18 and Llk = 42. Horizontal 
dashed line shows the point. Vertical dashed lines denote the 
jamming thresholds in the thermodynamic limit. 


located at a very well defined value in the 0-axes determining 
the jamming coverage for each k. The procedure in Fig. 1 was 
repeated for A:(= 3,4,5,6) (data do not shown here for reasons 
of space) and the results are collected in the second column of 
Table I. 

The high accuracy of the simulation data can be asserted by 
comparison with exact results for 9j presented in Ref. Jst]. In 
fact, for a fully aligned system, as studied here, the A:-mers are 
deposited only along an axis of the lattice, and the jamming 
problem reduces to the ID case. In this limit, 9j can be calcu¬ 
lated from Eq. O Si. The results are shown in the first col¬ 
umn of Table I. A remarkable agreement is obtained between 
theoretical and simulation data, validating the applicability the 
method introduced here to calculate jamming thresholds. 

In the next section, we will analyze the dependence of the 
jamming coverage on the size of the deposited objects for a 
system of straight rigid rods on a simple cubic lattice. 
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SIMULATION RESULTS: RANDOM SEQUENTIAL 
ADSORPTION OF STRAIGHT RIGID RODS ON A SIMPLE 
CUBIC LATTICE 

The procedure described in Section II was applied to cal¬ 
culate the jamming coverage corresponding to a system of 
straight rigid rods on a simple cubic lattice. A:-mers of length 
k = 2... 64 were considered, and the finite-size study was 
performed for lattices ranging between (36 x 36 x 36) and 
(1152 X 1152 X 1152) sites. In addition, n = 10^ runs were 
carried out for each pair k and L. These values of the param¬ 
eters allowed us to obtain very accurate measurements of the 
jamming thresholds, with errors less than 0.008% in all cases. 

In Fig. 2, the probabilities Wl( 0) are shown for two dif¬ 
ferent values of k {k = 2 and A; = 5 as indicated) and L/k{= 
18,24,30,36,42). From a simple inspection of the figure (and 
from data do not shown here for a sake of clarity) it is ob¬ 
served that; (a) for each k, the curves cross each other in a 
unique point Wj); (b) those points do not modify their numer¬ 
ical value for the different k used (ranged between k — 2 lo 
k = 64). In the cases of Fig. 2, = 0.509(9) for k — 2, and 

Wl = 0.505(9) for k = 5; (c) those points are located at very 
well defined values in the 0-axes determining the jamming 
threshold for each k, being Oj = 0.918388(16) for k = 2, and 
Oj = 0.736061 (12) for k = 5; and (d) dj decreases for increas¬ 
ing k-mer sizes. 

In the case of dimers (k = 2), the value of 9j obtained in Fig. 
2 differs significantly from the value 9j = 0.799(2) reported 
by Tarasevich and Cherkasova Q. Due to the methodology 
used in this contribution, our estimate of 9j is expected to 
be more accurate than that reported previously. In order to 
reinforce the validity of our result, typical configurations of 
dimers adsorbed on a simple cubic lattice are shown in Fig. 
3: (a) low-density state (0 w 0.4), (b) intermediate-density 
state (0 Ri 0.8), and (c) jamming state (0 Ri 0.918). Different 
colors represent different directions. Clearly, the value of 9j 
informed in Ref. ll^ corresponds to an intermediate state and 
is far from the limit density of the system. 

The procedure of Fig. 2 was repeated for k ranging between 
2 and 64, and the results are shown in Fig. 4 (solid circles) and 
second column of Table II. The points corresponding to k =32, 
48 and 64 were calculated for three relatively small values 
of L/k (12, 18, 24). The jamming coverage decreases upon 
increasing k. At the beginning, for small values of k, the curve 
rapidly decreases. However, it flatten out for larger values of 
k and Anally asymptotically converges toward a definite value 
as k ^ oo. 

Following a similar scheme to that proposed in Ref. 
for straight rigid rods on 2D square lattices, the curve of 9j 
vs. k was fitted to a function 

0,(k)= 0.4045(19) + ^:^-^:^ (k>6), (4) 

being 9j{°°) = 0.4045(19) the result for the limit coverage of 
a simple cubic lattice by infinitely long k-mers. 

Figure 4 includes also the behavior of the percolation 
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FIG. 4: The thresholds dj (solid circles) and Op (open circles) as a 
function of k. The solid line corresponds of the best fit of the jam¬ 
ming values as indicated in the text. Inset: Ratio Op/Oj as a function 
of k. 


threshold Op as a function of k for the system studied here 
(open circles). The corresponding numerical values, obtained 
from Ref. 1241] . are tabulated in the third column of Table 
II. The percolation threshold is a decreasing function of k. 
Combining the results of Oj and Op, the ratio Op/Oj was cal¬ 
culated (see inset in Fig. 4 and fourth column of Table II). 
As for 2D square lattices, the ratio between percolation and 
jamming thresholds exhibits a non-universal behavior. In this 
case. Op/Oj decreases in all range of k. 


CONCLUSIONS 

Random sequential adsorption of straight rigid A:-mers de¬ 
posited on a simple cubic lattice has been studied by numer¬ 
ical simulations and finite-size analysis. Several conclusions 
can be drawn from the present work. 

On the one hand, a new theoretical scheme to determine 
jamming thresholds was introduced here. The method relies 
on the definition of the probability Wl{0) that a lattice com¬ 
posed of Lx Lx L elements reaches a coverage 0. The value 
of Oj can be obtained from the crossing point of the curves 
of Wl for different lattice sizes. Comparisons with rigorous 
analytical results demonstrated the method’s accuracy. 


Table II: Values of Oj, Op and Op/Oj (as indicated in the text) 
for isotropic rods on a simple cubic lattice and k ranging 
from 2 to 64. 



Isotropic rods 

k 

Oj [Hi's] 

Op Ref. [24] 


2 

0.918388(16) 

0.2555(1) [23] 

0.2777 

3 

0.838860(14) 

0.2129(1) 

0.2538 

4 

0.780344(16) 

0.1800(1) 

0.2307 

5 

0.736061(12) 

0.1555(1) 

0.2113 

6 

0.701346(13) 

0.1364(1) 

0.1945 

7 

0.673355(14) 

0.1218(1) 

0.1809 

8 

0.650282(11) 

0.1089(1) 

0.1675 

9 

0.630901(12) 

0.0990(1) 

0.1569 

10 

0.614384(15) 

0.0901(1) 

0.1467 

11 

0.600130(10) 

0.0831(1) 

0.1385 

12 

0.587696(11) 

0.0772(1) 

0.1314 

13 

0.576780(14) 

0.0714(1) 

0.1238 

14 

0.567044(11) 

0.0661(1) 

0.1164 

15 

0.558360(13) 

0.0632(1) 

0.1128 

20 

0.525676(18) 

0.0478(1) 

0.0909 

24 

0.507750(18) 

0.0411(1) 

0.0807 

32 

0.483360(34) 

0.0299(1) 

0.0620 

48 

0.456071(32) 

0.0191(1) 

0.0417 

64 

0.440655(34) 

0.0143(1) 

0.0318 


On the other hand, the jamming coverage dependence on 
the k-mer size has been reported. In the case of dimers (k = 2), 
Oj = 0.918388(16) was obtained. This value differs signifi¬ 
cantly from the value Oj = 0.799(2) calculated by Tarasevich 
and Cherkasova iEbIi . Due to the methodology used in this 
contribution, our estimate of Oj is expected to be more accu¬ 
rate than that reported previously. 

On the basis of the behavior of Oj vs. k in the range 
2 < k < 64, and the best fit to this curve, the expression 
0j{k) = 0.4045 -F2.63/k —5.2/k^ was found, being 0j{°°) — 
0.4045(19) the result for the limit coverage of a simple cubic 
lattice by infinitely long k-mers. 

Finally, and based on previous findings from our group 
0241] . the possible relationship between percolation threshold 
and jamming coverage was examined. The results showed a 
non-universal behavior for the ratio Op/Oj, which decreases 
upon increasing k. 
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